Calculation of Stark-induced absorption on the 6s6p 3 Pi — 6s 2 1 Sq transition in Hg 
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We carry out relativistic many-body calculations of the Stark-induced absorption coefficient on the 
254- nm 6s6p 3 Pi(F = 1/2) — 6s 2 \So line of 199 Hg atom, the effect considered before by Lamoreaux 
and Fortson using a simple central field estimate [Phys. Rev. A 46, 7053 (1992)]. The Stark- 
induced admixing of states of opposite parity opens additional Ml and E2 transition channels. We 
find that the resulting Ml-El absorption dominates over E2-E1 absorption. The value of the E2-E1 
absorption coefficient depends strongly on the details of treatment of the correlation problem. As 
a result, our numerical values differ substantially from those of the earlier central field calculation. 
Reliable calculation of this effect can enable a useful experimental check on the optical technique 
being used to search for a permanent electric dipole moment of the 199 Hg atom. 

PACS numbers: 31.15.A-, 32.10.Dk, 78.20.Ci 



I. INTRODUCTION 

The F = 1/2 —> F = 1/2 electromagnetic transition 
between two atomic states of opposite parity has nec- 
essarily the electric-dipole (El) character. However, an 
application of the external E-field £ s breaks the spher- 
ical and mirror symmetries of the atomic Hamiltonian 
and opens all multipolar transition channels. To the 
lowest order in £ s , the transitions are determined by 
the Ml (magnetic-dipole) and E2 (electric-quadrupole) 
channels. These effects modify the absorption coefficient 
of the atomic sample, the corrective Ml and E2 terms 
scaling linearly with the electric field [l], Q • 

Lamoreaux and Fortson @ have focused on a specific 
setup, relevant for the search of thepermanent electric 
dipole moment (EDM) of Hg atom [3|, |4[ (non- vanishing 
EDM would violate P- and T-reversal symmetries and 
be a clear signature of new physics beyond the standard 
model of elementary particles). They considered excit- 
ing the 254-nm 6s 2 1 Sb — * 6s6p 3 Pi transition of 199 Hg 
atom. This isotope has the nuclear spin / = 1/2. A laser 
resolves the hyperfine structure of the 3 P\ level. It reso- 
nantly drives transitions from a given magnetic Mp i sub- 
level of the Fi = 1/2 ground state to the Ff = 1/2 level 
of the excited state. Then, for the F t = 1/2 -> F f = 1/2 
transitions, the relative change in the absorption coeffi- 
cient a may be parameterized as 
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Here £l is the polarization vector and is the direction 
of propagation of the laser wave. (Fi) is the expectation 
value of the total angular momentum, i.e., the nuclear 
polarization in the ground Ji = atomic state. 

In the current 199 Hg EDM experiment |4j, the 254-nm 
transition is used to monitor the nuclear spin direction 
and thereby detect EDM-induced shifts in nuclear spin 
precession, which will be linear in an external electric 
field. The Stark interference effect on the 254 nm ab- 



sorption given in Eq. ([T]) also is linear in electric field 
£ s and depends upon the nuclear spin direction. A reli- 
able calculation of this effect can enable a useful check on 
the EDM method when the effect is measured under the 
same experimental conditions as in the EDM search [5j. 



II. EXPRESSIONS FOR ABSORPTION 
COEFFICIENTS 



The goal of this paper is to compute the atomic- 
structure coefficients a mi and ciei ■ One may qualita- 
tively understand the appearance of Ml and E2 admix- 
tures in Eq. ([1]) as follows. The Stark-induced transition 
amplitude in a laser field is composed from terms lin- 
ear in the interactions with the external E-ficld, — D • £ s 
(D being the dipole operator), and the driving 2 fe -pole 
laser field. We may recouple the products of the two ten- 
sors (D and the laser EM multipolar interaction); the 
resulting compound operators have the multipolarities of 
|fc-l|,fc,fe + l. For the Fi = 1/2 -> Ff = 1/2 transition, 
k would be limited to 1 and 2. The additional constraint 
imposed by the parity selection rule yields the Ml and 
E2 multipolar couplings. 

We derived the expressions for clmi and a^i using the 
multipolar expansion of the plane EM wave and the first- 
order perturbation theory in the Stark field for the wave- 
function. We employ a geometry where the quantiza- 
tion axis z is chosen along the /c-vector of the linearly- 
polarized laser. The DC Stark-field and the laser polar- 
ization are aligned along the a:-axis, and the atom has 
a definite value of F y in the initial state. This particu- 
lar choice of geometry is convenient for working with the 
most general relativistic expressions for the multipolar 
transition operators [f| . By evaluating Eq. ([T]) of Ref . @ 
in this geometry, we identify the following expressions for 
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the structure factors 



a Mi 
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(3) 



Unless specified otherwise, here and below we use the 
atomic units, h — |e| = m e — 1, and the Gaussian units 
for EM equations. Here (rij Jj||r||n/ J/) is (within a fac- 
tor of —1 corresponding to the charge of the electron) 
the conventional reduced dipole matrix element for the 
6s 2 x Sq — ► 6s6p 3 Pi transition and fcz, = lul/c is the mag- 
nitude of the wavevector of the laser. The quantities 
R M1 and R E2 are sums over a complete set of interme- 
diate states; these sums arise due to the Stark-induced 
perturbation and involve the static El operator and mul- 
tipolar AC couplings to the driving laser field. Each of 
the sums, R M1 and R E2 , may be further split into two 

A'l 1 1 _Z?2 

sums Syj. , subscript i or / indicating which of the 
states, initial or final, is Stark-perturbed, 



R 
R 



Ml 



Ml 



L2 = -lV2Sf 2 (l°)-2 x ^-Sf 2 (2 e ) 

o V -L O 



(4) 
(5) 



The argument of the sums 5 1 , 



M1/E2, T7r 
// 



(J*') indicates the to- 



tal angular momentum J and the parity n of the inter- 
mediate states as fixed by selection rules. Explicitly, the 
reduced sums arc 
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We employ the relativistic formalism for the multipolar 
transition operators Q( M1 / E2 \ Specific single-particle 
reduced matrix elements computed using Dirac orbital 
parameterization of Ref. Q are 

<*IIQ (bj) IIj> - tellO/lh-) 

r J '{G t {r)G 3 {r) + F,(r)F, (r)}dr, (6) 



(i\\Q (mj) \\j) = 



J + l 

r J {G,(r)Fj(r) +Fi(r)G 3 -(r)}dr. (7) 



In both expressions we used the long-wavelength approx- 
imation, as auiL <§C 1. In these expressions, G(r) (F(r)) 



are the large (small) radial components of the Dirac bi- 
spinor, k are the relativistic angular quantum numbers, 
and Cj{f) are the normalized spherical harmonics. 



III. ATOMIC-STRUCTURE FORMALISM 

Mercury atom has two valence electrons outside a 
closed-shell core and we start our calculations with the 
so-called frozen core (V N ~ 2 ) Dirac-Hartree-Fock (DHF) 
approximation. In this approximation, the core orbitals 
are obtained self-consistently, while excited (valence) or- 
bitals are subsequently generated by solving the Dirac 
equation in the resulting potential of the core. Such or- 
bitals correspond to the Hg + valence orbitals. They are 
used as a basis for the standard configuration interaction 
(CI) technique for two valence electrons (see, e.g. 0]). 
We refer to this approximation as CI-DHF. Further sig- 
nificant improvement of the accuracy of the calculations 
is achieved when the standard CI technique is combined 
with many-body perturbation theory (MBPT) to include 
correlations of the valence electrons with the atomic core 
(CI+MBPT). 

The CI+MBPT formalism has been discussed in a 
number of papers (see, e.g., [1, H, G3). The effective op- 
erator (self-energy, S) arising from the core polarization 
may be split into a single-particle, Si, and a two-particle, 
E2, part acting in the model space. Qualitatively, a field 
of the valence electron induces an electric dipole of the 
polarizable core: Si describes an interaction of the va- 
lence electron with the self-induced core dipole, while £2 
describes its interaction with the core dipole induced by 
the other valence electron. We compute the self-energy 
correction in the second order of MBPT for the residual 
Coulomb interaction. Effects of higher orders will be also 
included in a semi-empirical fashion, discussed below. 

We use the Brillouin-Wigner flavor of MBPT to 
avoid the "intruder-state problem" , when the virtual core 
excitations inside £2 become resonant with the states of 
the valence subspace. Finally, we emphasize that our 
computations are ab initio relativistic and employ the 
Dirac equation and bi-spinors throughout the entire cal- 
culation. 

We use the second-order MBPT to calculate the self- 
energy operators Si and £2 via direct summation over 
a complete set of single-electron states. This set of basis 
states is constructed using the B-spline technique fi~ll |. 
We use 40 B-splines of order 9 in a cavity of 40 Bohr 
radius. The same basis of the single-electron states is 
also used in constructing the two-electron basis states for 
the CI calculations. We employ partial waves I = — 4 
and the 14 lowest states above the core in each partial 
wave (si/2, P1/2, P3/2, etc.) for the valence CI subspace 
and I — — 5 and 30 lowest states in each partial wave 
for internal summations inside the self-energy operator. 

Higher-order correlations are also included in £ in a 
way similar to Ref. [l2l |. The Si operator depends on the 
symmetry of the valence orbital. Therefore, we have a set 
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of different £1 operators for Si/ 2 , P1/2, P3/2, etc. states. 
An analysis of the spectra of Hg (see Table Q| shows 
that accurate treatment of Si is most important for s- 
electrons, because the ground 6s 2 state and other states 
with s electrons come close to the core and therefore core- 
valence correlations must be sizeable for them. In con- 
trast, the core-valence correlations are much smaller for 
more diffuse p and d orbitals. It turns out that the best 
accuracy is achieved if the all-order [13j operator is 
used for the s electrons, the second-order £1 is employed 
for the p electrons and no Si is included for d and higher 
waves. 

Higher-order contributions to £ 2 are included semi- 
cmpirically via screening factors which modify Coulomb 
integrals of the second-order £ 2 (see Ref. [l2| for details). 
The values of these factors are f = 0.9, fi = 0.72, / 2 = 
0.98, h = 1, f A = 1-02 and f 5 = 1.02. These values are 
found from comparing second-order and all-order £1. 

Finally, we further rescale the £1 operator for the s and 
p electrons to fit the experimental spectrum better. The 
rescaling coefficients are X s = 1.0961 and A p = 0.8675. 
We use the same A p for pi/ 2 and P3/2 waves. Note that 

X s > 1 because high-order effects, included in £1 for the 
s electrons, significantly reduce its value. On the other 
hand, A p < 1 because the second-order MBPT always 
overestimates the correlation correction. 

The resulting energies are listed and compared with 
experiment in Table [T] A typical deviation from the ex- 
perimental values is in the order of 100 cm" 1 . Even after 
the scaling, the disagreement remains, as the number of 
fitting parameters is limited. 

TABLE I: Experimental and theoretical energy levels of Hg 
(in cm -1 ). 



State J Exp. [14] Theory 



6s 2 










.000 


-13 


.79 


6s6p 


3po 





37645 


.080 


37458 


.26 






1 


39412 


.300 


39312 


.86 






2 


44042 


.977 


44265 


.45 


6s6p 


lpo 


1 


54068 


.781 


54180 


.72 


6s7s 




1 


62350 


.456 


62171 


.92 


6s7s 


L s 





63928 


.243 


63672 


.24 


6s7p 


3po 





69516 


.66 


69211 


.87 






1 


69661 


.89 


69385 


.18 






2 


71207 


.51 


70094 


.95 


6s7p 


lpo 


1 


71295 


.15 


71189 


.34 


QsQd 


L D 


2 


71333 


.182 


71295 


.01 


6s6d 


J 'D 


2 


71396 


.220 


71353 


.26 



The diagonalization of the CI+MBPT Hamiltonian 
provides us with the atomic wavefunctions and energies. 
While the wavefunctions already have correlation correc- 
tions built in, evaluating matrix elements requires addi- 
tional inclusion of the so-called screening effect. This 
effect arises already in the first order in the residual 
Coulomb interaction and describes a re-adjustment of the 
core orbitals in response to an externally applied field. 



We incorporate the screening in the framework of the all- 
order many-body technique, the random-phase approxi- 
mation (RPA). The RPA formalism (see, e.g., Ref. [11]) 
describes a linearized response to an oscillating perturba- 
tion. In this regard, while evaluating the reduced sums, 
we need to fix the driving RPA frequency for the entire 
set of matrix elements Q^ Ml ^> and Q^ E2 ^ at the photon 
frequency, u>l = Ef — Ei. However, for the dipole matrix 
elements (Stark mixing), the RPA frequency u> = 0. 

The evaluation of the sums S requires summing over 
a complete set of intermediate atomic states \n n J n ). 
We use two approaches: (i) direct summation over 
states (this implies explicit computation of the atomic 
states and evaluation of matrix elements), and (ii) 
the Dalgarno-Lewis method. In the Dalgarno-Lewis 
method [161 ] , the summation is reduced to solving the in- 
homogeneous Schrodinger (Dirac) equation (setup is sim- 
ilar to Ref. [13] )• As an illustration, consider evaluation 
of the sum Sj' n . It may be represented as 

sf'^MWQ^mf), 

where lumps corrections to the atomic wave func- 

tion of the final state due to the external field. It satisfies 
an inhomogeneous equation 

(H eS -E f )\6* f )=-r\n f J f ), (8) 

where if e ff is the effective CI+MBPT Hamiltonian of the 
atom. 



IV. NUMERICAL RESULTS 

As an illustration of the CI+MBPT methodology, we 
start with calculations of the El matrix element and en- 
ergy interval for the 6s 2 ^So — Gs6p 3 Pi transition. This 
matrix element normalizes the Stark-induced corrections 
to the absorption coefficient, Eq. ([T]). The theoretical re- 
sults at various levels of approximation and a comparison 
with the experimental values are presented in Table UTI 
We observe that the core-polarization (£1) has a substan- 
tial effect on the energy interval, leading to an improve- 
ment in the theory-experiment agreement. While the 
CI+£i value of the matrix element perfectly agrees with 
the experiment [l8|], such an agreement is fortuitous: in- 
cluding the screening correction to the Hamiltonian (£ 2 ) 
increases its value by a factor of 1.6; only the additional 
inclusion of the RPA screening and semi-empirical scaling 
moves the theoretical value into a 10% agreement with 
a 2%-accurate experiment. We find such an accuracy 
acceptable, as ab initio matrix elements of the intercom- 
bination (spin- forbidden) transitions are known [19j to be 
very sensitive to many-body corrections, the entire val- 
ues being accumulated due to the relativistic effects. On 
the other hand, the matrix elements of spin-allowed tran- 
sitions are stable with respect to inclusion of the MBPT 
effects (see, e.g., Ref. [Hj]). We will return to the evalu- 
ation of the accuracy of our calculations later. 
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TABLE II: Energy interval AE (in cm -1 ) and the reduced 
electric-dipole matrix element (R, a.u.) for the 6s 2 So — 
6s6p 3 Pi transition in Hg atom in various approximations. 



Approximation 


AE 


R 


CI-DHF 


31028 


0.405 


CI+Ei 


37441 


0.453 


CI+Ei+E 2 


37623 


0.716 


CI+E1+E2+RPA 




0.577 


as above but with all-order £2 


36947 


0.512 


as above but with scaled Ei (F: 


inal) 39313 


0.503 


Experiment [14, 181 


39412 


0.453(8) 



The Stark-induced correction to the absorption coef- 
ficient involves two channels, Ml and E2. We start by 
discussing the more involved oe2 calculations. We need 
to compute two sums, Sf 2 (J* = 1°) and Sf 2 (J£ = 2 e ). 
We carry out calculations (i) by direct summation over 
the 10 lowest-energy intermediate states of each symme- 
try (1° and 2 e ) and (ii) by using the Dalgarno-Lewis 
method. The latter method is equivalent to summing 
over infinitely many intermediate states. Both calcula- 
tions use the most sophisticated CI+MBPT approxima- 
tion (i.e., CI+E1+E2+RPA with semi-empirical scaling). 
The results are presented in Table 11111 An examination 
of contributions reveals that there are substantial can- 
cellations inside individual sums. This leads to an en- 
hanced sensitivity to correlations. For example, consider 
the value of the Sf 2 (J* — 1°) sum truncated at the 10 
lowest-energy levels. It changes from —35.6 to —28.23 
(Table Hill) wnile progressing from the CI+Ei+RPA 
to the full CI+MBPT treatment. Additional cancella- 
tions occur when the reduced sums are combined into 
the quantity R E2 = -\y/2Sf 2 (1°) - 2^Sf 2 (2 e ) a 
24.13 - 32.84 = -8.71. Notice that this value is sev- 
eral times smaller than the properly rescaled value of 
the largest contribution in Table HTT1 These cancellations 
may lead to a poor accuracy of our resulting absorption 
coefficient 

a E2 = -4.39 x 10 -3 a.u. = -0.0853 x 10 _8 /(kV/cm) . 

This result was obtained using the ab initio matrix ele- 
ment from Table III1 Notice that there is a phase ambigu- 
ity originating from atomic wavefunctions for sums S and 
the normalizing dipole matrix element. However, when 
these quantities are combined in Eq. (|2l3p . the ambiguous 
phase factors cancel out. In our particular computation, 
the sign of the dipole matrix element (6s 2 ^Sbl |r| |6s6p 3 Pi) 
is fixed by the first entry of Table IIII1 

We proceed to a comparison with results of Ref. 
These authors use a simplified approach in which a true 
many-electron problem is reduced to a set of single- 
electron problems. For the E2 interference they use the 
Dalgarno-Lewis summation method based on the DHF 
orbitals of the optically active valence electron. The 
LS coupling scheme was used in calculations. Their P- 
(D-)channel results correspond to our 1° (2 e ) values. 



Unfortunately, Ref. 0] contains a number of numeri- 
cal mistakes in calculations of a,E2 coefficient, hinder- 
ing a comparison. For example, for the P-channel, using 
Eq. (39) of Ref. and their numerical values we obtain, 
&E2,p = —0.96 x 10~ 8 cm/kV which is an order of mag- 
nitude smaller than the published value. Similarly for 
the D-channel, based on Eq. (46) and numerical values 
of Ref. @, we find aE2,D = 0.1 x 10~ 8 cm/kV, a factor 
of 20 smaller than the published value. We present a 
detailed comparison with (revised) values of Ref. for 
the two symmetries of intermediate states (1° and 2 e ) 
in Table HVl Although of the same order of magnitude, 
the individual contributions differ by signs. The most 
probable reason for the disagreement is the sensitivity to 
particulars of the treatment of correlations. For example, 
in computation of the D-channel contribution, Ref. [2j ne- 
glected intermediate states of the of 3 D symmetry. Hg is 
a heavy atom, and according to our table IIII1 omitting 
the triplet contributions would increase o.ei,d by a fac- 
tor of three. There is a remarkable cancellation between 
individual channels (Table IIV|) : our final result becomes 
an order of magnitude smaller than the recomputed value 
(a E2 = -0.87 x 10- 8 cm/kV) of Ref. 0. 

TABLE III: Breakdown of contributions to the reduced sums 
for the E2 Stark-induced transition. All quantities are in 
atomic units. 

sr W. = i°) 

n n J n (6s 21 g ||r||n n J n ) (n n J n \ \Q (E2) \\QsQp 3 Pi) contribution 



6s6p J Pi -0.503 7\949 22^29 

GsGp^i -2.956 -4.535 -54.41 

6s7p 3 Pi -0.037 -5.460 -0.63 

6s7p 1 Pi 0.674 -1.647 3.42 

6s8p 3 Pi -0.005 1.839 0.03 

GsSp^i -0.286 1.652 1.35 

6s9p 3 Pi -0.063 -2.875 -0.50 

6s9p x Pi 0.314 -0.269 0.23 

Sum(10) -28.23 

Dalgarno-Lewis, Sum(oo) —25.60 
Sf {JZ = 2 e ) 

n n J n (6s6p 3 Pi||r||n n J„) {n n J n \ |Q (E2) ||6s 2 ^jp) contribution 

6s6d L D 2 1.570 -6.963 75.18 

6s6d 3 L» 2 2.360 3.248 -52.61 

6s7d 1 D 2 0.576 -6.153 20.67 

6s7d 3 D 2 1.354 1.661 -13.10 

GsM^ 0.263 -2.857 4.10 

6s8d 3 L» 2 -1.363 0.422 3.14 

6s9d 1 D 2 -0.284 3.654 5.47 

6s9d 3 D 2 -1.281 -0.549 -3.71 

Sum(10) 39H 

Dalgarno-Lewis, Sum(oo) 44.97 



Fortunately, while ciE2 has a poor accuracy, it turns 
out to be much smaller than <Xmi , which, as shown below, 
can be computed reliably. There are two reduced sums 
to evaluate, Sf n (1°) and Sf 11 (l e ). Non-relativistically, 
the magnetic-dipole operator is diagonal in the radial 
quantum numbers. This means that the only substan- 
tial contributions arise in the sum (1°). Indeed, we 
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TABLE IV: Comparison of different multipolar contributions 
to the Stark- induced absorption coefficients cimi an d aE2 in 
l/(kV/cm). The first column gives the character of the multi- 
pole, the second column lists values of Ref. @], and the third 
column gives the results of our computation. The notation 
a[b] — a x 10 b is used. 



Contribution 



Ref. [21 



This work 



E2, 1° 
E2, 2 e 
E2, total 

Ml, l e 
Ml, 1° 
Ml, total 



-9.6[-9] a 
1.0[-9] a 
-8.7[-9] a 



7.8[-9] 
7.8[-9] 



2.4[-9] 
-3.2[-9] 
-0.85[-9] 

-0.13[-9] 
8.9[-9] 
8.86[-9] 



"Values recomputed by us based on data of Ref. [2j ; there 
are errors in numerical evaluations of Eq.(40), Eq.(46), and 
Eq.(47) of Ref. @]. See text for details. 



find from our fully relativistic analysis 



Sf L (l°) ra 0.0285, 
(l e ) « 0.0004. 



gMl 



The two dominant matrix elements entering 
Sf n (l°) are (6s6p 3 P 1 \\Q( M1 )\\6s6p 3 P 1 ) and 
(6s6p | IQ^ 1 ) | |6s6p ^i) . Both matrix elements 
may be estimated non-relativistically (e.g., one could 
use the Lande formula for the first matrix element). 
Further, the term involving the \6s6p 3 Pi) state is 
larger by roughly a factor of 5 than the contribution 
from the singlet state. As a result, the uncertainty 
in evaluating Sf 11 (1°) comes from the dipole matrix 



element entering this contribution, the already discussed 
(6s6p 3 Pi||r||6s 2 %). Incidentally, this is the very 
same matrix element that normalizes the absorption 
coefficient, so it cancels out in ajfi- Therefore, with 
about 25% accuracy 



o>M\ 



2( a P 1 HQ( M1 )j| 3 P 1 ) 



Ef - Ei 



1.19 x 10~ 8 /(kV/cm) . 



where we used the non-relativistic value 
( 3 -Pi||<3 (M1) || 3 -Pl) = (|) V&a. Our full-scale Dalgarno- 
Lewis relativistic CI+MBPT calculation results in 

OMi = 0.886 x 10" 8 /( k V/cm) , 

and is consistent with the non-relativistic estimate. From 
the preceding discussion, it is clear that our theoretical 
value is stable with respect to neglected many-body cor- 
rections. Ref. @ arrived at the result a^/i = 0.780 x 
10 _8 /(kV/cm). This differs by 12% from our estimates. 

Finally, we combine the contributions of the Ml and 
E2 interferences. We note that the poorly known E2 
contribution is fortunately suppressed by a factor of 10 
compared to the Ml coefficient. We find 

ami +a E 2= 0.80 x 10~ 8 /(kV/cm) . 
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